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Abstract 

The  dispersion  of  car  exhaustion  in  street  canyon  is  studied  by  LES  and  unsteady 
RANS.  Both  velocity  field  and  distribution  of  pollutant  concentration  are 
computed.  The  influence  of  the  geometry  of  building  blocks  on  the  pollution  is 
investigated. 

1.  Introduction 

The  prediction  of  the  pollutant  dispersion  is  of  great  significance  for  the  air 
quality  in  urban  area  ll'  2].  The  source  of  the  pollution  comes  from  the 
emission  released  by  the  motor  vehicles  at  the  street  floor.  The  large  eddy 
simulation  is  used  to  predict  the  unsteady  flow  as  well  as  the  dispersion  of 
pollutants  in  the  street  canyon.  Unsteady  RANS,  so-called  VLES  is  also 
applied  to  this  flow  and  the  results  show  that  LES  and  VLES  give  similar 
results  in  2D  cases  and  they  are  acceptable  for  the  environmental  problem. 


2.  The  governing  equation 

The  governing  equations  for  LES  can  be  written  as 


9u,  d(?qM  j)  _ 1 dp  d2uj  dxy 

dt  d xj  p d Xi  dxidxi  d x j 


(1) 

(2) 


3c+_  dc  _D  d2c  dMj 
dt  1 d Xj  dxidxl  d Xj 


(3) 


in  which  ty  and  M,-  are  the  subgrid  stress  and  mass  flux  respectively. 
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The  total  shear  stress,  i.e.  the  subgrid  stress  plus  molecular  viscous  stress,  is 
closed  by  the  RNG-based  Smargorinsky  model [3]  as  follows 

*,J  + Tff  = ~(veff  Sij  + /3)  (4) 

in  which  v eff  = v + v, . The  effective  viscosity  is  defined  as 

veff  =v[l  + H(y*veff/v3 -cT  (5) 

with  Sjj  ^i^Uj/dxj  +duj/dxi)/2  , vs  -ifmgVl^fJ2S~S~  and  H(x)=  x when 

x>0,  h(x)-0  otherwise.  The  V is  the  grid  volume  and  Crng  =0.157,  C=1.00. 

The  subgrid  mass  flux,  including  the  molecular  diffusion,  is  closed  by  a gradient 
model  such  that 

Mi=D,%L  (6) 

9*/ 

where  the  turbulent  diffusion  coefficient  is  defined  as  D,  =veff  /Sc,  , in  which 

Sc,  is  the  turbulent  Schmidt  number  and  assumed  to  be  1 .0. 

When  we  use  k-e  model  to  close  the  subgrid  Reynolds  stress  and  mass  flux,  it 
is  so-called  VLES.  The  k-  and  e-equations  with  standard  coefficients  are  ignored 
and  can  be  found  elsewhere  [4].  Same  formulae  as  eq.  (6)  is  used  for  average  mass 
flux. 

The  finite  volume  method  is  applied  for  the  numerical  solution  with  QUICK 
scheme  and  the  time  advance  is  second  order  Runge-Kutta  integration. 

3.  The  testing  cases 

3.1  Two  dimensional  street  canyon 

The  geometry  of  the  street  canyon  is  illustrated  in  Figure  l.The  upstream 
wind  speed  which  is  assumed  to  be  varied  with  time  period  of  8 minutes,  i.e. 

Uoo(y,t)=C/0(y)sin(2w/480).  The  pollutant  is  located  at  the  street  floor  and 

assumed  as  a constant,  presumably  the  pollutant  concentration  equals  to  1 .0  at  a 
local  spot  and  zero  outside  of  the  spot. 
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Figure  1 The  illustration  of  the  street  canyon,  • is  the  spot  of  the 

pollution 


The  computational  domain  and  grid  meshes  are  sketched  in  Figure  2 for 
B/H=  1.  The  mean  velocity  and  Reynolds  stress  profiles  are  given  at  the  inlet  and 
the  fully  developed  flow  condition  is  posed  at  the  outlet.  The  wall  function  is 
used  at  the  solid  walls  and  free  stream  condition  is  given  on  the  upper  boundary. 


Figure  2 The  computational  domain  and  the  grid  meshes 


3.2  Street  canyon  on  the  cross  road 

The  cross  road  is  illustrated  in  Figure  3.  A symmetrical  condition  is  added  in 
respect  to  the  plane  at  y=0.  The  car  exhaustion  is  simplified  as  a line  source  in 
the  middle  of  the  street  and  it  located  at  y axis  at  the  floor  of  street  canyon,  i.e. 
z=0.  The  geometry  and  flow  parameters  of  the  testing  cases  are  listed  in  Table  1. 

l i I i i 


> V 


Figure  3 The  plane  layout  of  street  canyon  on  cross  road 
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Table  I The  geometry  and  flow  parameters 


Case  A1 

Case  A2 

Case  B1 

Case  B2 

B/H 

1 

3 

1 

3 

Reynolds  number 

10b 

10b 

10b 

10b 

Turbulent  model 

LES 

LES 

k-e 

k-e 

3D  cross  road 

yes 

no 

no 

no 

H (meter) 

30 

30 

30 

30 

4.  The  results 
4.1  Case  A1  and  B1 

The  flow  fields  at  two  typical  time-steps  are  presented  in  Figure  4 and  5.  It 
clearly  shows  a vortex  in  the  street  canyon  that  transports  the  pollutant  outside  of 
the  canyon.  The  iso-contours  of  concentration  are  shown  in  Figure  6 and  7. 


(a) t=240s  (b) t=480s 

Figure  4 The  flow  field  for  B/H=  1 by  £-8  model 


(a)  240s  (b)  480s 

Figure  5 The  flow  field  for  B/H= 1 by  LES 


Note  that  the  velocity  and  concentration  fields  are  phase  averaged  on  six 
periods  in  LES.  The  flow  patterns  inside  the  street  canyon  are  nearly  same  in 
both  LES  and  VLES,  however  more  vortices  outside  of  the  canyon  in  LES, 
particularly  at  the  corns  of  the  building. 
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Results  of  the  concentration  distributions  are  in  good  agreement  between 
LES  and  VLES  inside  the  street  canyon  and  the  concentration  is  greater  on  the 
lee  side  of  the  front  building  and  smaller  at  the  front  of  the  back  building.  The 
concentration  distributions  are  different  greatly  outside  the  street  canyon  between 
LES  and  VLES,  this  is  due  to  the  different  flow  fields. 


(a)  240s  (b)  480s 

Figure  6 The  iso-contours  of  concentration  for  B/H=  1 by  model 


(a)  240s  (b)  480s 

Figure  7 The  iso-contours  of  concentration  for  B/H=  1 by  LES 


4.2  Case  A2  and  B2 

The  results  for  B/H=3  are  presented  in  Figure  8 and  Figure  9.  It  is  clearly 
shown  that  the  pollutants  are  transported  more  quickly  out  of  the  street  canyon 
for  the  wider  spacing  of  the  building  block. 

Compared  with  Figure  4 and  5,  Figure  8 and  9 clearly  show  that  the  bigger 
vortices  are  formed  in  the  wider  street  canyon.  It  is  useful  for  transporting  the 
pollutants  outside  of  the  canyon  and  reduced  the  concentration  of  pollution  on 
both  front  and  back  buildings,  as  shown  in  Figure  10  and  1 1 . 

Compared  between  LES  and  VLES  results  it  is  found  that  the  main  large 
vortices  are  nearly  same  inside  the  canyon  but  flow  patterns  are  different  in  the 
outside. 


(a)  240s  (b)  480s 

Figure  10  The  iso-contours  of  concentration  for  B/H= 3 by  model 


(a)  240s  (b)  480s 

Figure  11  The  iso-contours  of  concentration  for  B/H= 3 by  LES 
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4.3  Integral  flux 

Integral  flux,  shown  in  Figure  12,  is  a useful  design  parameter  in  building 
engineering  and  defined  as  the  integration  of  the  vertical  mass  flux. 


(a)  LES  (b)  k-e 

Figure  12  The  integration  flu x^w'c'dx 


Both  LES  and  k-e  model  show  the  same  periodicity  as  the  free  stream 
velocity.  On  average  the  agreement  between  LES  and  k-E  model  is  good  but  the 
peak  values  are  different.  It  is  due  to  the  different  flow  and  concentration  fields. 

4.3  Three  dimensional  dispersion  at  cross  road 

The  typical  velocity  and  concentration  fields  are  presented  in  Figure  1 3 and 
14.  The  velocity  fields  are  in  fairly  agreement  between  LES  and  k-e  models, 
however  the  difference  is  obvious  outside  of  the  street  canyon. 


(a)  LES  on  cross  section  at  y=100  (b)  k-e  on  cross  section  at  y=100 

Figure  13  Velocity  field  by  LES  and  k-e 
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The  similar  situation  appears  in  the  concentration  field  as  shown  in  Figure  14. 
Inside  the  street  canyon  the  concentration  distributions  are  quite  similar  but  big 
difference  occurs  outside  of  the  canyon. 


(a)  LES  at  y=100  (b)  k-B  at  y=100 

Figure  14  Iso-contours  of  concentration  in  3D  street  canyon  by  LES  and  k-B 

5.  Concluding  remarks 

Both  LES  and  VES  give  similar  results  for  flow  and  concentration  fields 
inside  the  street  canyon  in  both  2D  and  3D  cases.  However  LES  has  better 
resolution  in  flow  field  than  VLES  and  is  expected  to  have  better  prediction  for 
more  complicated  street  canyon.  On  the  other  hand  LES  spend  more  computation 
time  than  VLES,  5-8  times  more.  Therefore  the  computational  efficiency  of  LES 
must  be  improved  and  the  models  for  unsteady  RANS  should  be  also  consider  in 
order  to  satisfy  the  complex  3D  flows. 
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